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A fundamental problem of forced magnetic reconnection has been solved taking into account 
the plasmoid instability of thin reconnecting current sheets. In this problem, the reconnection is 
driven by a small amplitude boundary perturbation in a tearing-stable slab plasma equilibrium. It 
is shown that the evolution of the magnetic reconnection process depends on the external source 
perturbation and the microscopic plasma parameters. Small perturbations lead to a slow nonlinear 
Rutherford evolution, whereas larger perturbations can lead to either a stable Sweet-Parker-like 
phase or a plasmoid phase. An expression for the threshold perturbation amplitude required to 
trigger the plasmoid phase is derived, as well as an analytical expression for the reconnection rate in 
the plasmoid-dominated regime. Visco-resistive magnetohydrodynamic simulations complement the 
analytical calculations. The plasmoid formation plays a crucial role in allowing fast reconnection 
in a magnetohydrodynamical plasma, and the presented results suggest that it may occur and have 
profound consequences even if the plasma is tearing-stable. 
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I. INTRODUCTION 

In 1985, Hahm and Kulsrud published a seminal paper 
[1] on a basic problem of forced magnetic reconnection, 
the so-called “Taylor problem” (after the physicist J.B. 
Taylor who first proposed it). In this model of forced 
reconnection, a tearing-stable slab plasma equilibrium 
is subject to a small-amplitude boundary perturbation 
that drives magnetic reconnection. Adopting a resistive 
magnetohydrodynamic (MHD) description of the plasma, 
Hahm and Kulsrud showed that after an initial phase in 
which a current sheet builds up, the forced reconnection 
process evolves according to a linear resistive phase and 
then a nonlinear Rutherford regime [2] . 

Successively, in 1992, Wang and Bhattacharjee [3] 
showed that under certain conditions the reconnection 
process in the Taylor model passes into the nonlinear 
regime according to a Sweet-Parker-like evolution [3] in¬ 
stead of a Rutherford one. These two scenarios have been 
successively investigated for the case of a visco-resistive 
MHD plasma by Fitzpatrick In particular, he deter¬ 
mined an expression for the threshold perturbation am¬ 
plitude required to trigger a Sweet-Parker-like evolution. 

Since the pioneer work by Hahm and Kulsrud, the 
Taylor problem has not been free from controversies [Bi¬ 
lls] and has continued to be of great interest in plasma 
physics [ni434j . In fact, Taylor’s model has consider¬ 
able implications for laboratory plasmas (such as those in 
tokamaks, in which resonant magnetic perturbations can 
lead to magnetic reconnection) as well as astrophysical 
plasmas (such as those characterizing the solar corona, in 
which magnetic reconnection can be forced by the motion 
of the photospheric footpoints). However, despite the 


numerous works related to this problem, the nonlinear 
possible scenarios of the magnetohydrodynamical Taylor 
model have not changed from those found by Hahm and 
Kulsrud [Tj, Wang and Bhattacharjee [5], and Fitzpatrick 

la¬ 
in this article, extending our preliminary results [35) . 
we give evidence that a different nonlinear scenario is 
possible. We show that under certain conditions, after 
the initial phase in which a very elongated current sheet 
forms at the resonant surface, there is a sudden transi¬ 
tion to a fast reconnection regime characterized by the 
disruption of the current sheet due to the development of 
plasmoids. In spite of the fact that it is known from quite 
a long time that elongated reconnecting current sheets 
can become unstable to the formation of plasmoids |36F 
I39j , this possibility has not been taken into account in the 
theory of forced magnetic reconnection in Taylor’s model. 
In other contexts, instead, there has been a growing in¬ 
terest in the investigation of the link between plasmoid 
formation and fast reconnection [301 - 131] . 

Here, we extend previous analysis of the Taylor prob¬ 
lem considering also the plasmoid-unstable regime. We 
discuss all the scenarios of the collisional magnetohydro¬ 
dynamical Taylor model within a unified framework and 
we examine the conditions under which these scenarios 
occur. Both resistivity and viscosity are taken into ac¬ 
count, since both may play an important role in labo¬ 
ratory as well as in astrophysical environments [611163) . 
We perform both analytical and numerical calculations, 
showing that the development of plasmoids causes the on¬ 
set of a fast reconnection regime. Therefore, the present 
analysis is important to understand the evolution of the 
collisional reconnection process from the current sheet 
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formation to fast magnetic energy conversion rates. 


II. BASIC EQUATIONS AND TAYLOR’S 
PROBLEM 

We consider a plasma described by the two- 
dimensional (2D) reduced MHD model |M1- Within this 
description the magnetic and velocity fields can be ex¬ 
pressed in terms of flux and stream functions as 

y. , (1) 

^{x,y,t)=e^x\/(l){x,y,t), (2) 

where is the unit vector along the ignorable coordi¬ 
nate. The flux and stream functions are governed by the 
resistive Ohm’s law 

{dt+v-V)'il; = -'qj^+Eo, ( 3 ) 

and the vorticity equation 

(5t -I- V • V) = B • -I- . (4) 

Here, jz = — is the electric current density and ujz = 
is the plasma vorticity, both in the z direction. The 
kinematic viscosity can be recognized as v, whereas rj 
stands for the electrical resistivity. The term 
represents an externally applied electric field, which is 
required to maintain the equilibrium magnetic field in 
the presence of a small but finite electrical resistivity. 

This system of equations is dimensionless, with the 
lengths normalized to a convenient scale length L, the 
magnetic field to a convenient scale field strength Bq, 
and the time to the Alfven time ta = L/va, where va is 
the Alfven speed based on Bq. Therefore, the normalized 
resistivity corresponds to the inverse of the Lundquist 
number rj = Dn/{ vaL), being the magnetic diffusiv- 
ity, while the normalized kinematic viscosity is linked to 
the magnetic Prandtl number by means of the relation 

= P-mV- 

We consider a slab plasma which is periodic in the y 
direction and is bounded by perfectly conducting walls 
at a; = ±1 (note that the normalizing length L is chosen 
to be the half distance between the walls). Moreover, we 
consider a tearing-stable equilibrium given by 

i/'^°^(x) = —a;^/2 , = 0 . (5) 

This gives also that By^^\x) = x, jz^°^(x) = 1, and 
Bx^^Hx) = Vx^^^x) = Vy^°^(x) = ujz^^\x) = 0 (the nor¬ 
malizing field strength Bq is chosen to be the in-plane 
equilibrium magnetic field at \x\ = 1). 

Now we suppose that the conducting walls are subject 
to a sudden displacement 

Xuj = ±1 T 2(t) cos(/c?/) , (6) 


where k = 2 t: jhy and S(t) = Sq (l — — (t/r)e“*/’’’) 

for t > 0 [2^ [24] . The perturbation time scale r is as¬ 
sumed to be significantly larger than the Alfven time 
Ta (thus the plasma can be considered in magnetostatic 
equilibrium everywhere except near the resonant surface 
at X = 0) but much shorter than any characteristic re¬ 
connection time scale. Furthermore, the displacement is 
assumed to be small compared to the half distance be¬ 
tween the walls, i.e. Sq <C 1. Consequently, the bound¬ 
ary conditions at the flux conserving walls can be writ¬ 
ten as V’(±l,y,t) = -l/2-4'(f)cos(fcy), j^(±l,?/,t) = 1, 
(/)(±1, y, t) = ±d/dt{E{t))sm{ky)/k and a; 2 (±l, y, t) = 0. 
Note that 4'(f) = E!(t) since 
reason we have also that = Sq. 

In Taylor’s problem, as well as in tearing mode the¬ 
ory, it is convenient to divide the spatial domain into 
two regions: an “outer region”, where non-ideal effects 
and plasma inertia can be neglected, and a thin “inner 
region” centered on the resonant surface at x = 0, where 
resistivity, viscosity, and plasma inertia can be impor¬ 
tant. Hence, in the outer region the vorticity equation 
reduces to 


B-Vj2=0. (7) 

Since the boundary perturbation is relatively small, we 
can write V’out = where 'ip^^^{x, y, t) is the per¬ 
turbation in the form ^qJ;,j(x, y, t) = y t) cos{ky). 

Therefore, the linear approximation of Eq. 0 yields 

X (s^V’oul - = 0. (8) 

The solution of Eq. ^ which has parity 'ip^^{x,t) = 
''Ponti~x,t) and satisfies the boundary condition 
^oul(±l>^) = /dx\^^_^E{t) is [I] 


V’oul(a’>^) = V’oul(0,0 ( cosh lex - 


+4-(t) 


sinh \ kx 
sinh k 

(9) 

In the vicinity of the resonant surface, namely around 
X = 0, the previous expression can be described by a 
Taylor expansion as 


1 


^ontiXyt) = V’ouUO.i) + 2 ^V’out(i) |x| -f O (x^) 


( 10 ) 


with A'i/jout indicating the gradient discontinuity 

= A'oV’iul(0,i)+A;4'(t), 


AV’out(t) = dxpjiut 
where 


0+ 


— d 

'-'x Vout 


A' =_ 

° tanh k 


( 11 ) 


( 12 ) 


is the standard tearing stability parameter [65], which is 
negative in Taylor’s problem, while 


a: = 


2k 

sinh/c 


(13) 
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parametrizes the contribution of the external source per¬ 
turbation to the gradient discontinuity. 

In order to determine the explicit evolution of the sys¬ 
tem, we have to solve an initial value problem and to 
match the solutions obtained in the outer and inner re¬ 
gions. The early evolution of the reconnected magnetic 
flux can be calculated in the linear approximation, as 
was done by Hahm and Kulsrud for the case in which 
iz/r? < 1 [5 and by Fitzpatrick for the opposite case 
with ^ 1 [5]. In the next section we revisit their 
linear results, whereas in the subsequent sections we in¬ 
vestigate the nonlinear evolution by performing both an¬ 
alytical and numerical calculations. 


III. LINEAR THEORY 


The time dependence of the reconnected flux and the 
associated reconnection rate can be obtained by asymp¬ 
totically matching the inner and outer solutions in the 
overlapping domain. With this aim in mind, we now 
consider the vorticity equation and Ohm’s law in the 
inner region centered on the resonant surface cc = 0. 
Linearizing Eqs. a and 0 by considering perturba¬ 
tions in the form {x,y,t) = cos{ky) and 

4 >^^\x, y, t) = (j)^^\x, t) s\n{ky), and assuming that in the 
inner region <C 9^, we get 

+ kxcj)^^^ = , (14) 


= kxdl^^^^ + . (15) 


Performing the change of variable x = kx and then 
Laplace trasforming and Fourier transforming the per¬ 
turbed quantities as 


pOC 

' 4 ’{x,g)= / tjj^^\x,t) e~^*dt 

Jo 


and 


' 0 ( 6 ', 5 )= / 0(4,ff)e "^^dx, 


(16) 


(17) 


the set of equations (14)-(15) can be rewritten as 
yip + idgcj) = —rjk^O^tp , 


(18) 


gk^O^cj) = ik^dg{0'‘ip) — uk'^O'^fp. 


(19) 


Then, eliminating tp from Eq. (191 by means of Eq. (181, 
we obtain the layer equation 


d9 


0^ 

g + gk‘^ 0 '^ 


dej )) - g0^^ - vk‘^0^Jp = 0 . ( 20 ) 


In addition, the inner solution has to match the outer so¬ 
lution at (relatively) large a;, i.e. small 0. This boundary 
condition can be obtained by neglecting resistivity in the 
linearized form of Eq. ([^, and substituting ip^^}^{x,t) 
with its Taylor expansion around a: = 0. Then, after 
using the Fourier-Laplace transform, the boundary con¬ 
dition for the inner solution becomes 


0(^, 5) 11 '0out(0, 9) ^ -\-l + 0 (0)^ 

in the small 0 limit, where we have defined 


A(5)^ 


A0out(g) 

0out(O,5) 


( 22 ) 


(23) 


The quantity A.(g) can be evaluated by solving Eq. 


(20) subject to the boundary conditions (21) and (22). 


Then, A(g) may be used in the Laplace transform of Eq. 
(11) to obtain the “outer” magnetic flux 


0out(O,5) = 


AL4'(g) 

A(g) - A(, ■ 


(24) 


The reconnected flux can be calculated by considering 
the Laplace trasform of Eq (14), which yields 


gipix, g) = -x(p{x, g) + ■qk^dlip{x, g) . (25) 


By means of the Fourier transform we can rewrite this 
equation as 


pOO 

gip[x,g) = -x'(p{x,g)-gk^ i 0 '^ip{ 0 ,g)e"^^d 0 . (26) 


Then, using Eq. (18) to eliminate ip and evaluating this 


equation at a; = 0, we obtain the Laplace transformed 
reconnected flux 


0(0,g) = 


i0^ 


2r]k^ r 
9 Jo 9 + 


;de^d0. (27) 


Finally, the inverse Laplace transform of 0(0, g) gives 
us the time evolution of the reconnected flux 0*^^^(O,t), 
whose time derivative is a measure of the reconnection 
rate. 


In the following we solve the layer Eq. (20) and calcu¬ 


late the reconnection rate in two different linear regimes 
of the reconnection process, the inertial regime and the 
visco-resistive one. Note that here we consider a visco- 
resistive plasma which is characterized by plasma resis¬ 
tivity g that is smaller than kinematic viscosity v. The 
opposite case can be deduced following the same proce¬ 
dure. 


A. Inertial phase 


This equation must be solved subject to the condition 
lim0(6»,g) = O. (21) 

0—>-oo 


We first consider the inertial phase, which occurs for 
[5] t ^ Ti, with Ti = v~^l'^k~'^l^. In this regime, resistiv¬ 
ity and viscosity can be neglected in the layer equation, 
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which reduces to 


+ 



(28) 


The solution of this equation that satisfies the boundary 
conditions ( 21 ) and ( 22 ) is 


= fi/'out(0,5)^^ ■ (29) 


Substituting this solution into Eq. 
small 9 limit gives us 


( 22 ) and taking the 


A(5) = - 


kiT 

9 


(30) 


Furthermore, substituting the solution (291 into Eq. (27) 
and recalling that in this regime <C g, we get 


_ 277^/^_ 

^>(0,5) =-^V'out(0,5) ■ 


(31) 


Note that since \—2rjk‘^/g^\ 1, the inertial phase is 

a non-constant-i/) regime. We will see in the next sec¬ 
tions that this may have important implications in the 
nonlinear evolution of the reconnection process. 

We can now evaluate V'out(0j5) by means of Eq. (24). 


Neglecting the small initial transient r and also Aq with 
respect to A{g), we obtain 


V^(0,g) = -^A;vI;o. (32) 

TT 5 '^ 

Finally, the inverse Laplace transform of this expression 
gives us the time evolution of the reconnected flux 

tpiit) = = gk — 'i>of , (33) 

TT 

while the reconnection rate in the inertial regime can be 
evaluated as 


2A' 

= gk - -'Sot ■ (34) 

TT 

Note that in the inertial phase the current density 
( 0 , t) (which is proportional to the reconnection rate) 
increases linearly in time, whereas the width of the cur¬ 
rent channel <5 ^ g/fc [T] shrinks inversely with time since 
g . This leads to the formation of a narrow current 
sheet at the resonant surface. The thinning of this sheet 
towards singularity (for t ^ oo) is prevented by resis¬ 
tive diffusion, which becomes important as the current 
density piles up. 


B. Visco-resistive phase 

We now consider the visco-resistive regime, which oc¬ 
curs when [5] t ^ Tvr = and the magnetic 


island width w is less than the linear layer width (de¬ 
fined later in the paper). We can roughly follow Refs. 
[SilST] to solve Eq. ( [^ by asymptotic matching the 
solutions obtained for large 6 and relatively small 9. The 
large 9 domain represents an inner sublayer in real space, 
whereas the relatively small 9 region corresponds to an 
outermost sublayer in real space. 

We first consider large values of 9 such that 9 > 
(g/ry/c^)^/^. In this asymptotic interval, recalling that 
we are assuming v > g, Eq. (20) may be approximated 
and cast into the form 


dgcjyy — w = 0 , 


(35) 


where w ^ = gvk'^ and ( 0 , g) indicates the function 
that approximate ^{9,g) for 9 > {gjgk'^)^t'^. The ana¬ 
lytical solution of this equation can be easily found by 
transforming it into a familiar form. To this aim, let us 
introduce the change of variable (^> = 9^t‘^Xi which leads 
to 


^ d9^ 


9^- 

d9 


+ w-^9^]x = 0. 


(36) 


Then, with one more change of variable, precisely z = 
9^/ (3w ^), we can write the previous equation in the form 
of the modified Bessel equation 


i (Px . ^dx 

dz^ 




(37) 


with a = 1/6. The general solution of this equation is 
X = Cila{z) -\- C 2 Ka{z), where 1^ and are the a- 
order modified Bessel functions of the first and second 
kind, respectively, while ci and C 2 can be determined 
from boundary conditions. Since Ia{z) grows exponen¬ 
tially while Ka{z) decays exponentially, the solution of 
Eq. (35) that vanishes for 9 —)■ -|-c» is 




(38) 


In the opposite case of relatively small values of 9 such 


that 9 < (^i/gk^) 


4N-1/6 


d, 


Eq. (20) may be approximated as 
02 


g +gk^9^ 


de<j3< = 0 , 


(39) 


where (j)<i9,g) indicates the function that approximate 
(/)(0, g) in this asymptotic interval. The general solution 
of this equation is (/< = (- 5 / 0 - 1 - gk'^9)co + C 4 , therefore, 
by applying the boundary condition (221 we obtain 




0 


2'Kk 


(40) 

The quantity A ( 5 ) may be evaluated by matching 
the functions (/> and (/< in the overlapping interval 
{glgY^"^/k <9 < In this region becomes 


I -9- fn ^ f vk9A{g) 
=*2 ^out(0>3) ( — — -1 


(41) 
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while the expansion of (^> that match the previous solu¬ 
tion is 


'^> = * I V-out (0,5) (^6^/3 ^ - 1 j . (42) 

where F indicates the Gamma function. Therefore, the 
matching between </)> and d>< gives us 


A(5) =5' 


(43) 


where we have defined 
1 / 1/6 

T" ^ TT 




T (- 

and = 


r(s)' 


(44) 


A{g) allows us to evaluate the reconnected flux, which 
is equal to the outer flux in the visco-resistive regime . In 
fact, recalling that in this regime g rjk'^O'^, Eq. (27) 
gives us 


V'(0,5) = '0out(O,5) • 


(45) 


a non-constant-')/' behavior and the current sheet is sta¬ 
ble. This phase then gives way (on the time scale of 
resistive diffusion) to the Rutherford regime. 

iii) A “plasmoid phase” occurs if the magnetic island 
that enters into the nonlinear regime has a non-constant- 
■0 behavior and the current sheet is unstable. 

In the following we analyze all these nonlinear scenar¬ 
ios and show how they depend on both the microscopic {v 
and r]) and macroscopic (df, A(,, Aq and k) parameters of 
the Taylor problem. While the first two scenarios where 
discussed for the first time in the pioneer works by Hahm 
and Kulsrud [I] and Wang and Bhattacharjee [3] (and re¬ 
considered later by Fitzpatrick in the large magnetic 
Prandtl number limit), the plasmoid-dominated phase of 
the Taylor problem is entirely new. Since only part of the 
analysis of the visco-resistive plasmoid phase is strictly 
related to the Taylor paradigm, the applicability of the 
results is much more general. 


A. Rutherford phase 


Then, neglecting the small initial transient r, from Eq. 
(24) we get 


^^(0,5) = 


T* 5 ^ - ^o9 ’ 

which can be inverse Laplace transformed to obtain 

A(^, 

^0 


(46) 


^„,(t) = V'(i)(0,t) = ^ (e^o‘/”* - l) . (47) 


As a result of the fact that the forced reconnection 
process evolves initially according to the inertial phase 
(which is a non-constant-'0 regime), the system proceeds 
in the nonlinear regime in agreement to a Rutherford 
evolution only if the inertial phase gives way to a visco- 
resistive phase (which is a constant-i/; regime) before the 
transition to the nonlinear regime. Therefore, if we sup¬ 
pose that the transition from the non-constant-'0 to the 
constant-')/' magnetic island occurs at t ^ Tyr, the exis¬ 
tence of the visco-resistive phase requires that [S] 


Therefore, the reconnection rate in the visco-resistive 
regime is 


dt^vrit) = A), 



(48) 


IV. NONLINEAR THEORY 

The linear tearing mode analysis breaks down and non- 
linearities must be considered when the magnetic island 
width grows to a size comparable or exceeding the lin¬ 
ear layer width, which for a visco-resistive plasma corre¬ 
sponds to O nil IMl IB5] 

d^r, = 4:{'n-c^)^^^ pP ■ (49) 

In particular, we will see that three possible nonlinear 
scenarios may occur. 

i) A Rutherford evolution [2] takes place if the mag¬ 
netic island that enters into the nonlinear regime is char¬ 
acterized by a constant-)/) behavior. 

ii) A Sweet-Parker-like evolution [1] happens if the 
magnetic island that passes into the nonlinear regime has 


w{Tyr) < . 


(50) 


Recalling that the width of constant-')/' islands is linked 
to the reconnected flux by means of the relation 


w{t) = , 


(51) 


we can specify condition ( |50[ ) by evaluating the recon¬ 
nected flux at Tyy. Note that in Eq. (51) the contribu¬ 


tion of the equilibrium current density does not appear 
explicitly since in this problem, without loss of general¬ 
ity, we have = 1. An estimation of the 

reconnected flux at Tyr can be obtained from Eq. (471, 
which gives us 


(O' ~ ^ (e^o-”’-/- - l) . (52) 

Since Tyr A, the previous relation may be further 
semplified by using the Taylor series expansion of the 
exponential function. Hence, we can estimate 


w{Tyr) ^ 4 



(53) 
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Finally, condition (50) can be rewritten as 






(54) 


which tells us that the reconnection process proceeds in 
the nonlinear regime according to a Rutherford evolu¬ 
tion only if the magnetic perturbation amplitude is suf¬ 
ficiently small (i.e., for a sufficiently small displacement 


Eo = 'So/{d^^°>/dx\^^_^)). ^ 

We now proceed to determine the reconnection rate in 
the Rutherford phase, i.e., assuming that condition (54) 
is satisfied. In this case, since the constant-'?/; approxima¬ 
tion implies that ip{x,y,t) = -|- t) cos(fcy), 

the magnetic island geometry is entirely determined in 
terms of its width, whose evolution follows the Ruther¬ 
ford equation 


= +A;vl/(i)[iV (55) 

77 at \^/ 

where I = 0.823 [2J [55] • Here we may again neglect 
the perturbation time scale t since it is assumed to be 
much shorter than any characteristic reconnection time 
scale. Performing the change of variables t = and 
w = wjw;^^ with 


41 / A(4-o 

r7(-A')V(-A')’ 


(56) 


Wa = 



(57) 


It is now straightforward to evaluate the reconnection 
rate as 

- I for t « Ta , (62) 

(-^OI'^A V^a/ 

dtipRit) = — ^ t^T^. 

(63) 

Note that in the Rutherford regime the magnetic re¬ 
connection process is slow and depends strongly on the 
plasma resistivity. However, in the following we show 
that the evolution of the reconnection process may be 
very different for larger perturbation amplitudes. 


B. Sweet-Parker-like phase 

An elongated current sheet characterizes the initial 
stage of the nonlinear reconnection process if the transi¬ 
tion to the nonlinear regime occurs while the magnetic 
island has a non-constant- 7 /> behaviour. In fact, non- 
constant-?/; islands grow faster than the current can dif¬ 
fuse out of the reconnecting region, resulting in a strong 
current sheet at the resonant surface [Ij. This nonlinear 
scenario occurs if [5| 


w{Tyr) > S, 


I'T} ■) 


(64) 


namely if a constant- 7 /> magnetic island cannot form prior 
to the transition into the nonlinear regime. Following the 
arguments adopted to obtain condition (54), inequality 
(64) can be rewritten as 


we can rewrite Eq. (55) as 


dw 1 

di 


(58) 


From this equation we can see that w increases in time 
and reaches ih ^ 1 when t ^ Its integration gives 
that [Tj 


w = tanh(w -|- t) . 


(59) 


Then, the time evolution of the reconnected flux can be 
easily calculated in two interesting limits of Eq. (59). 
For t we may expand tanh(r() -|- t) into its Taylor 

series, which leads us to obtain 


^«(t)=(^y for t«.A. (60) 

In the opposite case t ^ Ta, we may neglect w with 
respect to t in the argument of the hyperbolic tangent, so 
that the time evolution of the reconnected flux becomes 


i’Rit) 


K'^0 

i-K) 


tanh^ 



for Tf^. 


(61) 


®"~ aSia = ®»' 


(65) 


which shows that the nonlinear evolution of the recon¬ 
nection process is characterized by an elongated current 
sheet as predicted in Ref. [3] if the perturbation ampli¬ 
tude is sufficiently large. 

We now analyze the evolution of the reconnection pro¬ 
cess when condition (651 is satisfied and also assuming a 
stable current sheet. An unstable current sheet is instead 
considered in the next subsection, together with the con¬ 
ditions for the instability. Here we calculate the recon¬ 
nection rate by performing a Sweet-Parker-like analysis 
in analogy with what was done by Wang and Bhattachar- 
jee 0, but taking into account also plasma viscosity ef¬ 
fects, which cannot be neglected if 1 / > ry. 

As is customary in this type of analysis, we consider a 
reconnecting current sheet (diffusion region) with length 
2Lcs and width 25 2Lcs- Then, assuming dy dx in 
the sheet and recalling that v • V?/; = 0 at the X-point, 
the reconnection rate may be evaluated from Eq. (|^ as 




( 1 ) 




= V 




( 1 ) 


X 




X 


B,, 


V- 


( 66 ) 
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where By^u indicates the magnetic field upstream of the 
current sheet. Since in this problem dz = 0, Faraday’s 
law V X E = — i9tB tells us that for a quasistationary 
magnetic field the ^ component of the electric field is 
uniform. Therefore, can also be evaluated from 

Ohm’s law in the ideal region just upstream of the recon¬ 
nection layer 




X 


^uBy^U • 


(67) 


Of course, indicates the plasma velocity upstream of 


the current sheet. Thus, from Eqs. (66) and (67) we can 
evaluate the width of the sheet as 





( 68 ) 


where we have estimated « Vd^jLcs from mass con¬ 
servation for an incompressible flow. 

We have now to evaluate the downstream velocity 
which may be done from energy balance considerations 
[H2] . Assuming a linear increase of the outflow along the 
current sheet and a Couette flow between x = S and the 
resonant surface x = 0, we obtain 


1 o 2 ^ 1 2 I I'd /2 

2By,u~ S'^ ' 


(69) 


Therefore, with the help of Eq. (68), we can find 


Vd ~ By ZL ( 1 -|- 

V 


- 1/2 


(70) 


In the limit ly/rj 1 we recover the standard Sweet- 
Parker result, namely that the downstream velocity cor¬ 
responds to the Alfven speed based on the reconnecting 
component of the magnetic field upstream of the diffu¬ 
sion region, Vd ~ By^y. However, for increasing values 
of viscosity the downstream velocity decreases due to en¬ 
ergy dissipation, and it approaches Vd « By^uirjjiyY^^ for 
vjr] > 1. 

We can now substitute the expressions for 5 and Vd 


into Eq. (66) to get 




X 


H3/2L-1/2 1/2 
^y,u ^cs '/ ' 7 y 


- 1/4 


(71) 


as found by Park et al. |62j in their extension of the 
Sweet-Parker model for visco-resistive MHD plasmas. 
However, here we want to go a step further by evaluat¬ 
ing By^y and Acs (supposed to be known in Sweet-Parker 
type analysis) as a function of the parameters of the Tay¬ 
lor problem. To this aim, we use Eq. (11) to express 


By^y as a function of the external forcing. Neglecting the 
small initial transient and the saturation term involving 
Aq, we can write By^y « A(,iE'o/2. Furthermore, numer¬ 
ical simulations indicate that the current sheet length 
is approximately one third of the periodicity length, i.e. 


2Lcb ~ By/3 = 27r/3fc. Therefore, from Eq. ( [tT] ) we 
obtain 


dt^IJsp = + 

X 3 \ Vj 


- 1/4 


(72) 

Finally, integrating this relation we obtain the time evo¬ 
lution of the reconnected flux 


1 / \ 

s == (i + () <. 

’ (73) 


Note that even if the reconnection rate in this modi¬ 
fied Sweet-Parker regime is higher than in the Rutherford 
one, it is still too slow to explain fast magnetic reconnec¬ 
tion events in plasmas characterized by very small resis¬ 
tivity and viscosity. We will see in the following that in 
such cases the plasmoid instability may resolve this issue. 


C. Plasmoid phase 


The analysis of the previous section relies on the as¬ 
sumption that the reconnecting current sheet remains 
stable during the evolution of the system. However, it 
should be pointed out that this is a strong assumption, 
which may be violated for small values of plasma resis¬ 
tivity and viscosity. In fact, despite the stabilizing effect 
of the outflow m, we might expect a tearing instability 
to arise in very thin reconnecting current sheets. A large 
number of numerical simulations have indeed confirmed 
this prediction [33 [371 SOI HIM EIHSS] ■ In particular, 
numerical simulations have shown that the tearing insta¬ 
bility arises if the aspect ratio of the reconnecting current 
sheet exceeds a certain threshold, i.e. when 


5 ^ Sq - ^cBc 5 


(74) 


being Cc the critical inverse aspect ratio, which was found 

to lie in the range [371 SSI SIl SSHSSl [ZHl ISP] 

This relation tells us that the analysis presented in the 
previous subsection applies if the width of the reconnect¬ 
ing current sheet is larger than a critical width, namely 
if (5 > Sc- Here instead we focus on the opposite case 
S < Sc- 


Let us first consider in more detail condition (74). As- 

1/2 

suming a critical current sheet of width Sc ~ {riLc/vd) ' , 
we can estimate Lc ~ tc'^rj/vd- Then, we may use Eq. 


(70) to get 


Lr 


e?H 


V,u 


1 + 


1/2 


(75) 


In our normalized units By^y = va , u , where va,u indicates 
the Alfven speed upstream of the diffusion region. There¬ 
fore, we can rewrite the condition for the tearing insta¬ 
bility of the (visco-resistive) reconnecting current sheet 
as 


Q _ BcsPA,u ^ q _ BcVAa 

O — > Oc — 



( 76 ) 




















FIG. 1: '^w = (Ag(red dashed line) and 5'c = 
C (fe/A(,) 77(1 + (blue solid line) as a function of the 

perturbation wave number k for rj = 10“®, u = 10“^ and 
C = 2(150)^. The critical wave number k* is approximately 
given by the intersection of the two curves. 


where Sc is the critical Lundquist number based on the 
critical length Lc- In the limit v ^rj\t follows that Sc ~ 
ISIIIB], whereas in the limit v ^ ij we obtain Sc « 
which agrees with the condition recently 
proposed by Loureiro et al. [54j. 

We are interested in a threshold condition for the 
perturbation amplitude rather than a condition on the 
Lundquist number based on the length of the current 
sheet. In order to achieve it, we recall that By^u ~ 
'Fo^s/ 2 and estimate Lcs k~^. Therefore, we obtain 
the condition 

'I^o>C—77(^1+ y ( 77 ) 


where C ~ 2e)r^ is a multiplicative constant. The frag¬ 
mentation of the reconnecting current sheet occurs when 
both conditions (651 and are satisfied. Hence, the 
plasmoid phase is triggered by the external perturbation 
when 


if ) 

Since 'Fry oc sinh/c and iFc oc sinhfc, there exists 

a critical wave number k* below which the reconnect¬ 
ing current sheet always leads to the plasmoid formation 
(see Fig. [^. Furthermore, since k* increases for de¬ 
creasing values of ry and v, the disruption of the current 
sheet is unavoidable in plasmas with very small resistiv¬ 
ity and viscosity. Note also that both and Tc de¬ 
crease with decreasing resistivity and viscosity, therefore 
the plasmoid phase may occur also for modest external 
perturbations if rj and v are particularly small. 

We can now evaluate the reconnection rate in the plas¬ 
moid phase as the rate of change of the flux reconnected 


at the main X-point. In fact, as nicely observed in Ref. 
[50] . only the open-flux parcels matter when evaluating 
the total reconnection rate, since a perfectly closed plas¬ 
moid does not carry any net reconnected flux. Therefore, 
summing and subtracting the contributions of all the re¬ 
connection layers one can find that the total reconnection 
rate is given only by the reconnection region at the main 
X-point. 

In a highly nonlinear regime the reconnection process 
is strongly time dependent, with plasmoids constantly 
being generated, ejected and merging with each other. 
However, we may assume a statistical steady state with 
a marginally stable current sheet located at the main X- 
point. In this case, the reconnection rate in statistical 
steady state can be evaluated as 


StV'p = 



(79) 


where (...) denotes time average and 



(80) 


Evaluating the velocity downstream of the marginal cur¬ 
rent sheet by means of Eq. (70), we can estimate 
Sc « r](l + iy/r/y'^^(ecBy^u)~^- Therefore, the resulting 
reconnection rate is 


S^t‘4^p 


I 1 + “ 


- 1/2 


(81) 


This expression shows us that when plasmoid-dominated 
reconnection develops in a visco-resistive plasma, the re¬ 
connection rate decreases with increasing values of mag¬ 
netic Prandtl number Pm = vlr\. This behavior is very 
different from what can be found when v rj. In this 
latter case, indeed, the resulting reconnection rate in sta¬ 


tistical steady-state becomes 


c R2 


and is 


therefore independent of the microscopic plasma parame¬ 
ters. However, as viscosity increases, this picture changes 
because the viscous energy dissipation leads to a decrease 
of the outflow velocity. 

Finally, in order to evaluate the reconnection rate only 
as a function of the parameters of the Taylor problem. 


we estimate By^u by means of Eq. (11). Therefore, we 
obtain the relation 


dti’p « ec(A(,«'o)" 


1 + - 
1 


- 1/2 


(82) 


which shows that the reconnection rate in the plasmoid- 
dominated regime depends strongly on the external forc¬ 
ing. A strong dependence was found also for the stable 
current sheet case, as can be seen from Eq. (72). This 


is because the external forcing increases the magnetic 
field upstream of the reconnecting current sheet. A sim¬ 
ilar observation was made also for laser-driven magnetic 
reconnection by Fox et al. [HU [HS], who explained the 
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fast reconnection rates observed in strongly laser-driven 
systems as the result of substantial upstream flux-pileup 
effects. 

Note that the disruption of the current sheet increases 
significantly the reconnection rate with respect to the hy¬ 
pothetical case of an elongated current sheet that would 
remain stable even for very small values of resistivity and 
viscosity. Indeed, Eq. (82) shows us that the reconnec¬ 
tion rate in the visco-resistive plasmoid phase depends 
on the magnetic Prandtl number = vjr] but not on 
the Lundquist number S, leading to the possibility of 
attaining fast magnetic reconnection even at very large 
Lundquist numbers. 

It is important to highlight that most of the presented 
theory applies to general reconnecting current sheets, ir¬ 
respective of being generated by a strong external forcing 
(large or a spontaneous process in a strongly un¬ 

stable equilibrium (large Ag), or a combination of both. 
Indeed, the condition ( [74| for the plasmoid instability of 
a reconnecting current sheet has general validity and ex¬ 
tends the previously proposed conditions [121021 [ 22 ] to 
arbitrary ratios vlrj. Similarly, also the relation (81) for 
the reconnection rate in the plasmoid-dominated regime 
is valid for both spontaneous and forced reconnection. In 
particular, this formula generalizes the previous results 
(i9tV’p « [ 22 ] for 1 / < ry and dtipp « ec.B 2 „(? 7 /iy)i /2 

[85| for y ^ r]) for arbitrary magnetic Prandtl number. 


V. NUMERICAL RESULTS 


We show now the different scenarios of the Taylor prob¬ 
lem by numerically solving the system of Eqs. (§ and 
(|^ with equilibrium and boundary conditions as speci¬ 
fied in Sec. |ll| To this aim, we employ a numerical code 
that splits all the fields into the time-independent equi¬ 
librium and an evolving perturbation, which is advanced 
in time according to a third order Adams-Bashforth algo¬ 
rithm. A compact finite difference algorithm is adopted 
to compute the spatial operations in the x direction (on 
a non-equispaced grid), while a pseudospectral method 
is used for the periodic direction. We adopt a space dis¬ 
cretization of 512x8192 grid points, so that the minimum 
step size in the x direction is dx^^^ = 1.6 x 10 “^, whereas 
in the y direction we have dy = 6.14 x 10“^. 

In the following, we consider the results of simulations 
with Ly = 16 tt (i.e., k = 1 / 8 ), r = 2 and different values 
of d/g, rj, V. The choice of the parameters is such to 
fall within the three scenarios presented in the previous 
section. Indeed, for a fixed perturbation wave number fc, 
the possible evolutions of the system can be organized in 
a three-dimensional “phase diagram” with 'I'g, 77 and v 
on the three axis (a four-dimensional “phase diagram” is 
needed to retain also the k dependence). Such a diagram 
can be constructed from the conditions derived in Sec. 
IV and summarized in Sec. VI Figs. [^a) and[^b) show 
two slices of this diagram (the Hahm-Kulsrud scenario, 
which occurs if dtg ^ 'i'w, is supposed to hold until 


(a) 



n 

(b) 



V 


FIG. 2: (Color online) Two dimensional slices of a “phase dia¬ 
gram” for forced magnetic reconnection in the magnetohydro- 
dynamical Taylor model. Fixed parameters are (a) k = 1 / 8 , 
v/r) = 10 , and (b) k = 1 / 8 , rj = 10 ~^. Three possible evolu¬ 
tions can occur: (i) Hahm-Kulsrud scenario mu, (ii) Wang- 
Bhattacharjee scenario OE], (iii) New scenario characterized 
by a plasmoid phase. These different scenarios of the forced 
reconnection process are summarized in Sec. jvll 


'I'o = dliy/3), which can be used to graphically organize 
the knowledge of the reconnection dynamics in a system 
with given plasma and forcing parameters. 

The typical evolution of the reconnection rate when 
the perturbation amplitude satisfies condition (54) (but 
is not so small as to prevent the system from entering into 
the nonlinear regime) is shown in Fig. In this case, 
a current sheet builds up at the resonant surface during 
the inertial phase, which lasts until t < Ti ~ 147. Then, 
after the magnetic island attains a constant-T/; behavior at 
t ^ Tyr ^ 684, the reconnection process evolves through a 
visco-resistive phase as long as the island width becomes 
comparable to the linear layer width at t ^ 7000. After 
this time the system proceeds nonlinearly according to 
a slow Rutherford evolution as predicted by Hahm and 
Kulsrud [T]. 
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FIG. 3: Reconnection rate as a function of time 

from a numerical simulation with 77 = 2 x 10“®, u = 2x 10~®, 
3'o = 4 X 10"^. 



FIG. 4: Reconnection rate as a function of time 

from a numerical simulation with 77 = 8 x 10~®, u = 8x 10~®, 
5'o = 4 X 10“^. 



0 100 200 300 400 500 

t 


FIG. 5: Reconnection rate as a function of time 

from a numerical simulation with 77 = 4 x 10“®, 7/ = 4 x 10~^, 
VI/q = 4 X 10“^. 



FIG. 6: (Color online) From the numerical simulation shown 
in Fig. contour plots of the out-of-plane current density 
jz(x,y) with the in-plane component of some magnetic field 
lines (white lines) superimposed at (a) t = 320 and (b) t = 
450. 


A different evolution is shown in Fig. In this 

case, the perturbation amplitude satisfies condition (65) 
but not condition ([77]). After the initial inertial phase 
{t < TiK, 93), the system passes into the nonlinear regime 
at t ^ 90 while the magnetic island is still characterized 
by a non-constant-i/) behaviour. As a consequence, the re¬ 
connection process evolves according to a Sweet-Parker- 
like regime and only subsequently a Rutherford regime 
takes place, as predicted by Wang and Bhattacharjee [3] 
A completely different scenario occurs when the per- 
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FIG. 7: Time-averaged reconnection rate in the plasmoid- 
dominated regime dttpp as a function of the magnetic Prandtl 
number Pm for 17 = 4 x 10“®. The theoretical slope Pm^^^ is 
shown for reference. 


turbation amplitude satisfies condition (781, as shown in 
Fig. § In this case the inertial phase lasts until the 
non-constant-i/i island width becomes comparable to the 
linear layer width, at t ^ 150 (well before « 543). 
Then the system proceeds into the nonlinear regime ac¬ 
cording to a Sweet-Parker-like evolution characterized by 
an elongated current sheet with slowly varying aspect 
ratio. The elongated current sheet that characterizes 
this stage of the reconnection process is shown in Fig. 
I^a). However, as the reconnecting current sheet grad¬ 
ually thins, the plasmoid instability sets in at t « 360. 
A sudden increase of the reconnection rate occurs when 
the first nonlinear plasmoids reach the width of the cur¬ 
rent sheet, which begins to break up at t « 410. Then, 
after t « 440, when the plasmoid chain is well developed 
(see Fig. [^b)), the reconnection rate seems to fluctuate 
around a constant mean value, as expected on the basis 
of previous numerical simulations of plasmoid-dominated 
reconnection [13 Sana EH E3 EH- The growth of the 
plasmoids is so fast that they start to merge before being 
ejected in the downstream region, and the simulation is 
stopped at t = 490, when the adopted truncated Fourier 
expansion can no longer resolve the large gradients in the 
y direction that are due to the coalescence of the plas¬ 
moids in this stage of the reconnection process. 

Although in our simulations we are not able to achieve 
a long statistical steady-state of the plasmoid phase, it is 
interesting to try to look at the dependence on the mag¬ 
netic Prandtl number of the reconnection rate in this 
phase. Performing numerical simulations by varying Pm, 
we obtain that the time-averaged reconnection rate in the 
plasmoid phase scales roughly as Pm when the mag¬ 
netic Prandtl number is large (see Fig. [^, in agreement 
with the theoretical scaling (82) in the limit Pm 3> 1. 
The reconnection rate for the simulation with Pm = 10 
deviates more from the theoretical predictions mainly be¬ 
cause of the numerical dissipation that adds up to the 


physical dissipation and that is more significant when 
the physical viscosity is lower. 


VI. SUMMARY 

In this work we have examined a fundamental prob¬ 
lem of forced magnetic reconnection, the so called Taylor 
problem, in which magnetic reconnection is driven by 
a small amplitude boundary perturbation in a tearing- 
stable slab plasma equilibrium. In particular, we have 
shown how the evolution of the reconnection process de¬ 
pends on both the external source perturbation and the 
microscopic plasma parameters. 

Three possible evolutions can occur in a magnetohy- 
drodynamical plasma. 

i) Hahm-Kulsrud scenario mm- if the perturbation 

is such that Iko ^ (^77)^^®, the reconnec¬ 

tion evolves initially through a linear inertial phase, then 
through a linear visco-resistive phase, and finally (if the 
perturbation is not excessively small) it proceeds in the 
nonlinear regime according to a Rutherford evolution. 

ii) Wang-Bhattacharjee scenario B [5]: if the per¬ 
turbation is such that (A(,fc^/^)“^(i/77)^^® ^ 'kg < 
C (fc/Ag) 77(1 + V , after the linear inertial phase the 
system passes into the nonlinear regime according to a 
Sweet-Parker-like evolution, which then gives way to the 
Rutherford regime. 

iii) Our scenario: if the perturbation is such that 'kg > 

(A(,fc^/^)“^(i/77)^^® and T'g > C (fc/A(,) 77(1 -|- , af¬ 

ter the linear inertial phase and an initial nonlinear evo¬ 
lution with a gradually evolving current sheet, the re¬ 
connection suddenly enters into a fast plasmoid phase. 
Note that the realization of a Sweet-Parker sheet before 
the onset of the plasmoid phase is not mandatory, since 
in systems with extremely large Lundquist numbers the 
plasmoid instability may be triggered prior to the attain¬ 
ment of such current sheet [sniisiiiHS] 

The plasmoid formation plays a key role in allow¬ 
ing fast magnetic reconnection by making the reconnec¬ 
tion rate independent of the Lundquist number, which is 
known to be extremely large in many astrophysical and 
laboratory plasmas. While for the case of negligible vis¬ 
cosity the reconnection rate is also independent of the 
magnetic Prandtl number, we have shown that when vis¬ 
cosity is not negligible « ec(A(,'kg)^(l -|- 77/77)”^^^. 
Therefore, the reconnection rate in the visco-resistive 
plasmoid phase decreases as the magnetic Prandtl num¬ 
ber increases. The reconnection rate depends also on the 
perturbation amplitude, whose effect is to increase the 
magnetic field upstream of the reconnection region. 

It is important to highlight that there exists a critical 
perturbation wave number k* below which the evolution 
of the system always leads to the plasmoid phase. Since 
k* increases for decreasing values of the plasma resistiv¬ 
ity and viscosity, the development of stable elongated cur¬ 
rent sheets is not attainable in many situations of interest 
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in astrophysical and laboratory fusion plasmas. There¬ 
fore, also modest perturbation amplitudes may lead to 
plasmoid-dominated reconnection in plasmas character¬ 
ized by very small resistivity and viscosity. 
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